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Abstract 

This paper pays attention to finding the free surface 
elevation during liquid sloshing at low excitation frequency. 
Experiments have been conducted at low frequency ratios to 
find the free surface elevation and the obtained results are 
reported here. The liquid behavior is captured by means of a 
video camera and the free surface oscillations of liquid at 
different excitation frequencies calculated by the image 
processing technique. The liquid height is observed by the 
intensity variation of the interface between the air and liquid. 
The location of maximum gradient has been considered as 
interface. The effects of excitation amplitude, frequency and 
liquid fill level have been studied. 
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Introduction 

The knowledge of natural frequencies for the liquid 
free surface is important in the design of liquid tanks 
subjected to different types of excitation. The dynamic 
behavior of a free liquid surface depends on the type 
of excitation and its frequency content. Excitation with 
frequencies in the vicinity of the lowest natural 
frequencies of the liquid motion is of primary practical 
interest. Civil engineers and seismologists have been 
studying liquid sloshing effects on large dams, oil 
tanks and elevated water towers under ground motion. 
Sloshing phenomenon is one of the major concerns in 
the design of liquid storage tanks undergoing ground 
excitation (Maleki and Ziyaeifar, 2008). Experimental 
studies can provide very useful and reliable results for 
some cases. However, they are generally very 
expensive. Moreover, there is a lack of high quality 
experimental data in the area of liquid sloshing 
(Akyildiz and Unal, 2005). Sloshing can be broadly 
classified into two types, namely self-induced sloshing 



and externally induced sloshing, based on how 
sloshing is generated inside its container. The main 
difference is that, in self-induced sloshing, the 
container is immobile whereas, in externally induced 
sloshing, the container will move owing to external 
disturbances. Several investigations on self-induced 
sloshing (Hu et al. 1999, Okamoto et al. 1991, Saga et al. 
2000, Singh et al. 2006, Nasar et al. 2009) in a 
rectangular tank with circulating flow have been 
conducted experimentally and numerically in the past 
20 years (Takizawa et al. 1992, Okamoto et al. 1993, 
Fukaya et al. 1996). Most of the studies on sloshing 
were limited to self-induced excitations. In the case of 
externally induced sloshing, an external force such as 
an earthquake or a sudden brake on a moving 
container vehicle acts as the excitation source for 
sloshing. The effect of sloshing in a moving liquid 
container causes safety and comfort issues such as 
undesirable forces on brake systems in liquid carrying 
tankers like LNG (Graczyk et al. 2006, Lee et al. 2007), 
high impact loads upon the containment system and 
structural parts. Further, under earthquake conditions, 
it creates fatigue, safety problems in nuclear reactors 
and power generating plants (Kimura et al. 1995). 
Detailed review of literatures can be found in Ibrahim 
(2005) and in Eswaran and Saha (2011). 

Experimental Plan 

In this investigation, the liquid behavior in an 
externally induced tank is captured by means of a 
video camera and the free surface oscillations of liquid 
at different excitation frequencies calculated by image 
processing technique. Finally, the statistical analysis is 
performed to show the probability density and power 
spectral density of signal. 
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Experimental Facilities and Procedures 

Figure 1 shows the schematic diagram of the 
experimental set-up, details of which have already 
been reported by Eswaran et al. (2009) and Panigrahy 
et al. (2009). The partially filled liquid tank of size 600 
mm x 600 mm x 600 mm is attached to a shaking table 
(915 mm x 930 mm) which can be moved to and fro by 
a crank arrangement and is driven by a DC motor 
(rated 1440 RPM, 220 V, 1 HP). As the required 
frequency of oscillation is too low, therefore, a 
rectifier-cum-controller is used to obtain the desired 
frequency of oscillation. The controller that converts 
the alternating current to direct current and feeds it to 
the DC motor with appropriate regulation is 
connected with the DC motor. It contains a full wave 
rectifier, which converts the 220V AC-1(|) supply to 
220V DC. 




Concrete Basement 



i 1 220 V, 1 Phase, 

All Dimensions are in mm AC Supply 



FIG. 1 BASIC EXPERIMENTAL SET-UP 

The crank arrangement is used in the set-up to 
transform the rotary motion to reciprocating motion 
which connects the DC motor and the platform of the 
shaking table. The water is used as a working fluid in 
the experiments. A video camera (Sony DCR-SR 300 
Handy cam Camcorder) is screwed onto an adjustable 
beam which is attached to the shaking platform with 
its lens focusing on the liquid free surface at the left 
corner of the tank as shown in Figure 2. The adjustable 
beam is moved up or down until a suitable slot is 
reached so that the camera focuses only on the liquid 
during its entire motion. 
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The camera is held at a distance of 300 mm away from 
the liquid free surface and focuses on the water free 
surface. Initially, the relation between the pixel and 
millimeter has been found by still images. The camera 
is then ready to capture the images of the air-water 
interface as shown in Figure 3. After adding a coloring 
agent (navy blue) in the liquid, an excellent contrast 
between the liquid and the white marker particles is 
obtained that gives a superior image quality with 
lesser noise. The motor is allowed to rotate whose 
speed can be varied using a speed controller. The 
speed of the motor is kept sufficiently low so that the 
sloshing occurs without much turbulence. The tank 
can be filled to any desired level and a scale attached 
to the side wall of the tank gives the height. The 
experiments have been conducted by varying the 
frequencies and the hJL ratio of the tank for different 
case studies. 




Binary im 



Image height 




FIG.2 CAMERA ARRANGEMENT IN THE SET-UP FOR 
FINDING THE FREE SURFACE ELEVATION 



Image Width 

FIG.3 CONVERSION OF RGB IMAGE INTO BINARY 
IMAGE 

Interface Location Technique for Capturing Free 
Surface Elevation 

The luminance contrast between the liquid and the air 
which indicates the location of the interface is 
determined by maximum intensity gradient method. 
The location of maximum gradient is considered as 
interface as described by Law et al. (1999) and Eswaran 
et al (2011). The videos captured by the video camera 
are then transferred to the computer, and further 
image analysis is carried out. The liquid height is 
observed by the intensity variation of the interface 
between the air and liquid. The intensity gradient at 
each point has been calculated. The location of 
maximum gradient has been considered as interface. 
The movement of interface can be analyzed by cross- 
correlation between these matrices. The free surface 
oscillations are plotted for different excitation 
frequencies. Figure 4 shows the typical pixel intensity 
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variation along a vertical line cutting across the 
interface region. The movement of the interface can be 
tracked by analyzing its displacement between two 
successive images. 
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FIG.4 INTENSITY VALUES FOR BINARY IMAGE 

All images are then available in the form of a vertical 
location matrix. The movement of interface can be 
analyzed by cross-correlation between these matrices. 
The standard deviation of extracted vertical locations 
denotes the interface fluctuation with time and the 
mean value gives the mean fluctuation for the 
corresponding excitation frequency and hs/L ratio. 

Tests to Ensure the Accuracy of 
Experimental Results 

Depending upon the type of liquid, nature of external 
disturbances on the tank, tank shape and hs/L ratio, the 
free liquid surface can experience different motions 
like simple planar, non-planar, rotational, irregular 
beating, symmetric, asymmetric, quasi-periodic and 
chaotic (Ibrahim, 2005). The amplitude of slosh, in 
general, depends upon the nature, amplitude and 
frequency of the tank motion, liquid-fill depth, liquid 
properties and tank geometry (Pal et al. 2002, Akyildiz 
and Unal, 2006). As the rectangular tank oscillates, 
different sloshing waves are created depending on the 
liquid depth and frequency of oscillations. In this 
study, the problem was restricted to liquid sloshing in 
a rectangular tank under surge (movement of tank 
only in horizontal direction) oscillations. The dynamic 
tank velocity (Vt) can be determined from the 
following equation, 



V t =-A t cocos(cot) 



(1) 



where At is the tank excitation amplitude and CO is the 
oscillation frequency. The amplitude of the tank 
excitation is 100 mm. For a given rectangular 
prismatic tank, the natural frequencies of the fluid 
depending on the fill depth are given by, 



(2) 



where L and d are the tank width and water depth 
respectively and n is the mode number (Akyildiz and 
Unal, 2005). In Table 1, different combination of liquid 
hs/L ratio, excitation frequency and excitation 
amplitude are shown for a total of thirty six cases. It is 
to be noted that at higher frequency, the turbulence 
level becomes higher. If the frequency of the tank 
motion is close to the natural frequency of the liquid, 
large sloshing amplitude occurs. Since the liquid 
elevation is calculated through image processing, the 
more turbulence of liquid should be avoided which 
will induce noises while processing. Thus, all the 
experiments have been conducted at low frequencies 
only. 

Uncertainty Analysis 

Uncertainty is a quantification of the doubt about the 
measured result. The maximum possible uncertainty is 
estimated from the minimum value of quantity 
measured and accuracy of the measuring system. The 
maximum error calculated while measuring the 
excitation frequency is 0.014 cycle/sec. The 
experiments are repeated and are found consistent. 
The detailed error analysis is discussed in Appendix A. 

Repeatability of Experiment 

In the experimental plan as mentioned in Table 1, the 
set 5 was used to the check repeatability of 
experiments. The repeatability of the slosh 
phenomenon was investigated by running the same 
experiment several times. The test is conducted by 
plotting the height versus time for case numbers from 
37 to 41 which is sketched in Figure 5. It is seen that 
the measurements are very close to each other. To get 
a measure of how close the experiments are, the 
standard deviation is calculated for each time instant. 
The calculation of the standard deviation a(t) is given 
by 



m = i\Y i (at)-rn(t)Y 



mi 



(3) 



(4) 



where ^. is the measured surface elevation of 
experiment i = 1,2 . . . ,5 and m(t) is the average of all 
experiments. The standard deviation &(t) is shown in 
Figure 6. The standard deviation is found smaller as 
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compared to the measured value, thus showing the 
repeatability of the tests. 
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FIG.5 REPEATABILITY OF TEST (CASE-37 to 41) 
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FIG. 6 STANDARD DEVIATION OF EXPERIMENTS (CASE-37 
to 41) 

Validation with. Wave Height Gauge 

Wave elevation can also be measured using wave 
height gauge. As shown in Figures 7 and 8, three 
distinct marks viz., upper, lower and midpoint are 
made on the probe wire for convenience. Initially, the 
lower mark is set to immerse in the tank liquid and the 
liquid level is adjusted to coincide with the middle 
mark. The difference in permittivity between water 
and air is used to measure the surface elevation in a 
moving tank. Then the gauge results are compared 
with the results obtained from interface location 
technique as shown in Figure 9. When the surface is 
oscillating, a thin layer of liquid is formed on the 
package wall. This thin layer of liquid moves down 
much slower than the surface of the liquid. However, 
the result matches well with allowable accuracy. 
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FIG.7 WAVE PROBE 





FIG.8 PERSPECTIVE VIEW OF WAVE PROBE AND 
ANALOGG TO DIGITAL CONVERTER 
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Interface location technique 

-a— Wave height gauge 




FIG.9 TANK 600 x 600 x 600 mm, FL=200 mm, FR=0.6, 
EXCITATION AMPLITUDE=60 mm 

Comparison of Experimental and Numerical Results 

In this section, the experimental results for two 
different frequency ratios are compared with the 
numerical simulation results. Numerical simulation 
procedures and results have already been reported by 
Eswaran and Saha (2009, 2010). A fully non-linear 
model for idealized 3-D waves in a numerical wave 
tank has been developed. A sigma ( <7 ) transformation 
is used to map the liquid domain onto a rectangle, 
such that the moving free surface in the physical plane 
becomes a fixed line in the computational mapped 
domain. The fourth order central difference scheme 
and the Gauss-Seidel point successive over-relaxation 
iterative procedure are used to capture the free surface 
wave profiles and free surface elevation plots of the 
fluid domain. 

The free surface elevation for both simulated and 
experimental results is compared during the motion of 
tank and a minor discrepancy has been observed. It 
can be observed from Figures 10 and 11 that both 
experimental and numerical results indicates the 
presence of a single traveling wave. Two cases are 
considered for such a validation by choosing two 
different frequency ratios. In the both cases, the 
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excitation amplitude and the fill ratio (fc/L=0.5) 
remains constant. 



Numerical (Sigma-Transformation) 

— — =— — Experimental (Interface Location) 




FIG. 10 FREE SURFACE ELEVATION FOR CASE-19 
(FR=27.4 %; hs/L =0.5; EA=30 mm) 

The excitation amplitude (EA) value has been chosen 
as 30 mm since better results can be expected with low 
excitation amplitude for the present experimental and 
numerical works. When the results are closely 
examined, one can observe that the numerical results 
clearly over scoop the experimental results for the case 
of high frequency ratio (Figure 12). But when the 
frequency ratio is low, the crest and trough points in 
both results are matched well. 
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FIG.ll FREE SURFACE ELEVATION FOR CASE-20 
(FR=36.6 %; hs/L =0.5; EA=30 mm) 

Results and Discussion of Free Surface 
Elevation 

Effect of Excitation Amplitude 

As the rectangular tank oscillates, different sloshing 
waves are created depending on the liquid depth and 
frequency of oscillations. These waves are asymmetric 



and at large amplitude the tank excitations are 
combined with the travelling waves. Liquid heights 
for different frequency ratios under low amplitude 
excitation (EA=30 mm) for cases 11 and 12 as per Table 
1 are shown in Figures 12 and 13 respectively. For the 
case of low amplitude excitation, the liquid height 
reaches the steady level quickly. It is observed that the 
slosh amplitude follows a sinusoidal pattern which 
means the initial increase of slosh amplitude with the 
increase in applied frequency and thereafter it 
decreases and again with the increase in external 
frequency the slosh amplitude increases. This is 
because of the sinusoidal nature of applied external 
excitation. 




40 60 

Time (sec) 

FIG.12 FREE SURFACE ELEVATION FOR CASE-11 
(FR=39 %; hs/L =0.35; EA=30 mm) 

As the amplitude of excitation is increased (EA=100 
mm), the liquid responds violently related to the fluid 
motion that occurs due to the occurrence of turbulence, 
hydraulic jump, wave breaking and 3-D effects. Figure 
14 shows the liquid height for different frequencies at 
100 mm excitation amplitude. It has been observed 
that the wave amplitude grows initially with the time 
up to a certain time period. After reaching this period, 
the wave amplitude remains constant and enters the 
fully developed region. In other words, in the first 
phase, the wave amplitude rises slowly from mean 
level up to certain period whereas in the second phase, 
a uniform oscillation can be observed. Similar trend is 
noticed for rest of the cases (case-28 to case-36) in set 4 
experiments of Table 1. Free surface oscillation or the 
instantaneous wave height (h) for first phase (i.e., to 
20 seconds) and for second phase (i.e., 100 to 120 
seconds) at tank left wall for case-28 is shown in 
Figures 15 (a) and (b) respectively, and their 
corresponding wave phase plane diagrams are 
depicted in Figure 16 (a) and (b). In each case, the 
instantaneous wave height is plotted for first phase 
and second phase time periods. 
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FIG.13 FREE SURFACE ELEVATION FOR CASE-12 
(FR=49 %; hJL =0.35; EA=30 mm) 
Phase I Phase II 
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FIG.14 FREE SURFACE ELEVATION FOR CASE-28 
(FR=29 %; ML =0.35; EA=100 mm) 
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FIG.15 FREE SURFACE ELEVATION FOR CASE-28 
(a) Phase I; (b) Phase II. 
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FIG.16 PHASE PLANE DIAGRAM FOR CASE-28 
(a) Phase I; (b) Phase II. 

Effect of Excitation Frequency and Fill Level 

As the tank moves, it supplies energy to sustain the 
liquid motion. The four possible types of waves are 
standing wave, traveling wave, hydraulic jump and a 
combination of these. Figures 17 and 19 show the 
effect of excitation frequency and effect of fill level. As 
the frequency increases, the standing wave transforms 
into a train of traveling waves of very short length, 
and the hydraulic jump occurs near first mode natural 
frequency. Further, with increase of frequency, the 
standing wave transforms into a train of traveling 
waves of very short length. Hydraulic jump, usually, 
takes place due to a small disturbance and appears 
over a range of frequencies near the resonant 
frequency (Akyildiz and Unal, 2005). With further 
increase in frequency, the jump will pass into a 
solitary wave. For deeper liquid, sloshing near the 
resonance is characterized by the formation of large 
amplitude standing waves. These waves are 
asymmetric in nature, and at large amplitude, tank 
excitations may be combined with the traveling waves. 
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TABLE 1 DETAILS OF THE EXPERIMENTAL PARAMETERS 



Set No. 


Case 


Excitation 
Amplitude 
(EA) 


Fill Level 
(FL) 


Excitation 
Frequency ( 00, ) 


First Mode 

Natural 
Frequency 
(&>!> 


FR 

(coJco^xlOO 






mm 


mm 


hJL 
ratio 


cycle/sec 


rad/sec 


rad/sec 


% 




1 


30 


120 


0.2 


0.3 


1.885 


5.34835 


35.24 




2 


30 


120 


0.2 


0.4 


2.513 


5.34835 


46.98 




3 


30 


120 


0.2 


0.5 


3.142 


5.34835 


58.74 


Setl 


4 


40 


120 


0.2 


0.3 


1.885 


5.34835 


35.24 


5 


40 


120 


0.2 


0.4 


2.513 


5.34835 


46.98 




6 


40 


120 


0.2 


0.5 


3.142 


5.34835 


58.74 




7 


60 


120 


0.2 


0.3 


1.885 


5.34835 


35.24 




8 


60 


120 


0.2 


0.4 


2.513 


5.34835 


46.98 




9 


60 


120 


0.2 


0.5 


3.142 


5.34835 


58.74 




10 


30 


210 


0.35 


0.3 


1.885 


6.41167 


29.40 




11 


30 


210 


0.35 


0.4 


2.513 


6.41167 


39.19 




12 


30 


210 


0.35 


0.5 


3.142 


6.41167 


49.00 


Set 2 


13 


40 


210 


0.35 


0.3 


1.885 


6.41167 


29.40 


14 


40 


210 


0.35 


0.4 


2.513 


6.41167 


39.19 




15 


40 


210 


0.35 


0.5 


3.142 


6.41167 


49.00 




16 


60 


210 


0.35 


0.3 


1.885 


6.41167 


29.40 




17 


60 


210 


0.35 


0.4 


2.513 


6.41167 


39.19 




18 


60 


210 


0.35 


0.5 


3.142 


6.41167 


49.00 




19 


30 


300 


0.5 


0.3 


1.885 


6.86364 


27.46 




20 


30 


300 


0.5 


0.4 


2.513 


6.86364 


36.61 




21 


30 


300 


0.5 


0.5 


3.142 


6.86364 


45.77 


Set 3 


22 


40 


300 


0.5 


0.3 


1.885 


6.86364 


27.46 


23 


40 


300 


0.5 


0.4 


2.513 


6.86364 


36.61 




24 


40 


300 


0.5 


0.5 


3.142 


6.86364 


45.77 




25 


60 


3UU 


0.3 


0.3 


1 QQC 
l.OOO 


6.86364 


27.46 




AO 


DU 


300 


0.5 


0.4 


2.513 




3A £1 




2/ 


6U 


300 


0.5 


0.5 


3.142 


6.86364 


A C r 7 r 7 
40./ / 




28 


100 


210 


0.35 


0.3 


1.885 


6.41167 


29.39 




29 


100 


210 


0.35 


0.35 


2.199 


6.41167 


34.29 




30 


100 


210 


0.35 


0.4 


2.513 


6.41167 


39.19 




31 


100 


210 


0.35 


0.5 


3.142 


6.41167 


49.00 


Set 4 


32 


100 


120 


0.2 


0.3 


1.885 


5.34835 


35.21 




33 


100 


120 


0.2 


04 


2.513 


5.34835 


47.26 




34 


100 


180 


0.3 


0.3 


1.885 


6.15000 


30.59 




35 


100 


180 


0.3 


0.4 


2.513 


6.15000 


40.81 




36 


100 


240 


0.4 


0.3 


1.885 


6.60861 


28.52 


Repeatability test 




37 


30 


300 


0.5 


0.8 


5.0265 


6.86364 


73.23 




38 


30 


300 


0.5 


0.8 


5.0265 


6.86364 


73.23 


Set 5 


39 


30 


300 


0.5 


0.8 


5.0265 


6.86364 


73.23 




40 


30 


300 


0.5 


0.8 


5.0265 


6.86364 


73.23 




41 


30 


300 


0.5 


0.8 


5.0265 


6.86364 


73.23 



(H = Height of the test tank; hs = Fill water level; L= Length of the tank) 
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FIG.17 EFFECT OF FREQUENCY FOR CASES - 13, 14 and 15 
(hs/L =0.35; EA=40 mm) 
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FIG. 18 EFFECT OF FILL LEVEL FOR CASES - 8, 17 and 27 
(FR « 49; hs/L =0.35; EA=30 mm). 

Figure 19 shows the bar chart representation of 
maximum and minimum wave heights for case 1 
through 4. The maximum and minimum wave height 
for first phase and second phase is represented as 
Imax, Imin and Smax, Smin respectively. It shows that 
the maximum wave height increases and minimum 
wave height decreases while increasing the frequency. 



ID Smax □ SMin 



200 
160 
120 
80 
40 




FIG. 19 BAR CHART REPRESENTATION OF MAXIMUM AND 
MINIMUM SLOSH DISPLACEMENT WITH FREQUENCY FOR 
CASES - 28 to 31(fc/L =0.35; EA=100 mm). (Fl - 1.885 rad/sec, F2 
- 2.199 rad/sec, F3 - 2.513 rad/sec, F4 - 3.142 rad/sec). 




Spectra of Wave Elevation 

The spectra of a wave elevation are computed by the 
Fast Fourier Transform (FFT) as shown in Figure 20, 
where S is the wave elevation spectrum. Figure 20 
shows the spectra of wave elevations at the tank left 
corner. It causes peaks in the power spectra and more 
peaks appear at different frequencies. The sloshing 
waves have been analyzed with the probability 
density and energy spectrum. 
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FIG.20 POWER SPECTRA OF WAVE ELEVATION AT THE TANK 
LEFT CORNER: (a) CASE-28; (b) CASE-29; 

(c) CASE-30; and (d) CASE-31 



Conclusions 

This paper is devoted to presenting the experimental 
procedure for finding the free surface elevation in a 
moving rectangular tank. The liquid free surface 
behavior at low frequency oscillation has been 
elaborately studied from experiments. Initially, the 
history of wave or long term behavior of free surface 
oscillation (0 to 120 seconds) is observed. Thereafter, to 
study in detail, two time phase periods are focused 
and instantaneous wave heights are plotted for each 
case. It has been observed that the wave amplitude 
grows initially with the time up to a certain time 
period. After reaching this period, the wave amplitude 
remains constant and enters the fully developed 
region. It is found that the maximum wave height 
increases with an increasing frequency, while the 
minimum wave height decreases with increasing 
frequency. It also causes the peaks in the power 
spectra to become larger and encourages more peaks 
to appear at different frequencies. The sloshing waves 
have been analyzed with the probability density and 
energy spectrum. It is clear that the trend of deviation 
of wave height increases with an increase in frequency. 
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This suggests that a stronger nonlinearity gives rise to 
the increase of the deviation. It is interesting to note 
that such an increase seems to indicate a fairly linear 
trend with frequency. Finally, the experimental data 
have been compared with the results obtained from 
the numerical simulation. 
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APPENDIX - Error Analysis 

Measured Quantities 

TABLE Al MAXIMUM POSSIBLE UNCERTAINTY FOR 
MEASURED QUANTITIES 



SI 
No 


Measur 
ement 


Device 


Accuracy 


Minimum 
Quantity 


Maximum 
Uncertainty 


1 


Time 


Stop Watch 


±0.5s 


120 s 


± 0.5 % 


2 


Pressure 


Piezo- 
resistive type 
transducer 


± O.Olbar 


5 bar 


± 0.2 % 


3 


Wave 
height 


Capacitance 
wave gauge 


± 0.2 mm 


0.05m 


± 0.4% 



In this experimental investigation the parameters time, 
pressure and liquid height are measured quantities. 
Stop watch is used to measure time and piezo-resistive 
sensor is used for pressure analysis. During motion of 
the tank, the liquid height is measured by interface 
location technique and traditional wave height gauge. 

Estimated Quantities 

The frequency parameter is estimated from the 
number of cycles completed per second. Here, the 

uncertainty is calculated as — = 0.42% using the 
procedure described by Kline and McClintock (1953). 
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